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We consider the K° — K° and B° — B° mixings in the MSSM with the two-Higgs- 
doublet scalar sector featuring explicit CP violation, and the Yukawa sector of type II. 
In the case of strong mixing between CP-odd and CP-even states the existence of light 
charged Higgs is allowed in the model. The mass splitting Am LS and the amount of 
indirect CP violation e are calculated. In the limit of effective low-energy approximation 
the nonstandard effects are shown to be negligibly small in Amis and e for the K°- 
mesons, being almost independent on the charged Higgs mass. However, for the — B^ 
and B a s — B® systems the effects of nonstandard physics are shown to be larger, limiting 
the MSSM parameter space. 
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1. Introduction 

Minimal supersymmetric extension (MSSM) of the standard model (SM) pQ in- 
volves many sources of CP violation besides the well-known CKM mixing. In the 
general case a number of complex parameters can be included in the MSSM leading to 
the observable effects of CP violation dependent on their phases [2j [3] in addition to 
the only CKM phase of the SM. 

The number of phases in the soft SUSY breaking terms can be reduced significantly 
using the assumption of trilinear parameters A t> b universality in the Higgs boson-squark 
sector and neglecting the gaugino phases. Restricted CP violation [I] of the Higgs 
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potential can be described by the only one universal phase ip = arg(^A b ) = a.xg(fiA t ) 
in addition to the CKM phase (/x is the Higgs superfield mixing parameter) and leads 
to a relatively simple picture of radatively induced CP violation in the MSSM two- 
doublet Higgs potential. A calculation of radiative corrections by means of the effective 
potential method leads to the complex parameters A 5 6i 7 depending on the products 
{iA t b with phases different by a factor of two, argA 5 = 2 argA 6 7 . Direct experimental 
reconstruction of phases could be performed in the processes of superpartners and 
Higgs bosons production at the colliders [5] while an indirect evidence could be given 
by measurements of electric dipole moments [6], mixings of the neutral mesons and 
meson decays [TJ. 

The MSSM two-doublet Higgs sector contains three neutral scalars and two charged. 
Radiatively induced complex parameters of the Higgs potential, when appear, lead to 
explicit CP violation and the three neutral scalars hi, h 2 , h 3 which are mixings [8], [9] 
of CP-even states h, H and CP-odd state A known in the limit of CP conserving two- 
doublet potential. The lightest scalar hi can have a mass substantially smaller than the 
LEP2 direct limit for the SM m# >114 GeV pU] being not observed at LEP2 energies 
because of the ZZh\ coupling suppression [11] by parameters of the h\, h 2 , h 3 mixing 
matrix which does not require specifically large values of the tg^—v^/vx parameter. 
Charged scalars H could be relatively light, with masses 50-100 GeV. In the framework 
of such MSSM scenario an additional box diagrams with charged Higgs exchanges 
could contribute noticeably to mass splittings and mixing parameters in the neutral 
meson systems. In this connection it is interesting to analyse MSSM implications 
to the neutral meson mixings, the only systems where effects of CP violation are 
experimentally observed. 

In the following first we consider the neutral K°- meson system where CP violation 
was discovered [12] . The approximations used in this study are also applied to the 
Bg, B° D mesons. In section II we briefly describe the MSSM Higgs sector and the 
Yukawa sector of type II (model MSSM II) and approximations used in the SM for 
calculation of the meson mass splittings Amis and the CP violation parameter e of 
the neutral i^-meson system. In sections 3,4,5 we calculate Amis an d £ including 
the one-loop box diagrams with one or two charged Higgs boson exchanges using the 
low-energy approximation. Some consequences are discussed in section 6. 
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2. Model MSSM II 



Yukawa sector. The Yukawa sector Lagrangian in a model with n-doublet ($ n ) 
structure of the scalar sector can be written in the form 

plus hermitian conjugated term. Here gf„ n , gf™ — are the Yukawa couplings of up 
and down quarks with scalar doublets which are 3 x 3-matrices in the flavor space with 
matrix elements wich are, generally speaking, complex, {u' iL , d' iL } are the left doublets 
of up and down quarks, u'j R are right singlets of up quarks and d'j R are right singlets 
of down quarks (i, j =1,2,3 define the quark generation number). 

The two-doublet Yukawa sector of type II [13] does not include the gfj 2 , gfj 1 terms, 
so $i doublet generates only the up quark masses and <3> 2 doublet only the down quark 
masses. The MSSM II Yukawa sector Lagrangian is given by 

-1% = g^Q'i^mU + g^Q'iL^d'jR + M. sec. + h.c. (2) 

In terms of physical fields (diagonalized Lagrangian) the neutral scalars A, H, h 
participate only in the flavor- conserving quark interactions, the charged scalars take 
part only in the flavor-changing transitions. The flavor- changing neutral currents do 
not appear. 

Scalar sector. The general hermitian renormalizable SU(2) ®C/(l)-invariant MSSM 
potential at the m top energy scale has the form 

+A 1 ($ t 1 $ 1 ) 2 + A 2 (^$ 2 ) 2 + A 3 ($ t 1 $i)(<l» t 2 $ 2 ) + A 4 ($ t 1 $ 2 )($ t 2 <l»i) + 
+ ^($t$ 2 )($t$ 2 ) + ^($t$ l) ($t$ l) + 

+A 6 ($ t 1 $ 1 )($ t 1 <l> 2 )+ A 6 ($l$i)(^$i) + A 7 (^$ 2 )($ t 1 $ 2 )+ A 7 2 )(^$i). 

At the Msusy mass scale (energy of the order of sparticle masses) Ai 7 parameters 
are real and expressed through the electroweak SU(2) ® U(l) couplings g\ and ^ 2 [H]. 
At the energy scale m top large radiative corrections to the scalar self-interactions can 
appear from the box and triangle one-loop diagrams with scalar quark exchanges, 
generating complex As,.. 7 parameters. In the restricted complex MSSM mentioned in 
the Introduction, these corrections depend on four parameters, the trilinear universal 
A tt b, Higgs superfield parameter /1, SUSY mass scale Msusy and the universal phase 
<P = arg(^A t)6 ). 
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The diagonalization procedure for the general two-doublet potential, leading to 
Higgs boson mass eigenstates and their self-interactions in the physical basis, can be 
found in [H [15]. Three of the eight scalar states from two complex doublets (the 
Goldstone modes G° and G^) are eaten by W ± - Z°- bosons, giving two charged 
scalars H and three neutral scalars hi, h2, h% without definite CP-parity which are 
the mixtures of h, H, A. The charged Higgs boson mass is given by [H [T5] 



m 2 H ± 



2 , 2 

m w + m A 



-(ReAA 5 - AA 4 ), 



(4) 



where my/ is the W boson mass, m A is the CP-odd scalar mass defined at the zero 
phase, v = 2mw/g2, and the effective parameters AA^s at the one-loop are [I] 
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The Yukawa couplings h t = ^^k, h b = and the X tb stands for 
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Radiative corrections to the effective potential defined by ([5]) and at the scale 
[f| are especially large at moderate Msusy (several hundreds of GeV) and large 
/i, A tjb parameters (of the order of 1 TeV). The scalar quark sector is strongly cou- 
pled. For the case of substantial phase arg(yuA t b ) the ^1,2,3 mass splittings at mod- 
erate H ± mass of 150-180 GeV are as large as 10-50 GeV. Unacceptable regions 
of the MSSM parameter space exist, where close to zero or not posi- 

tively defined. Note that the correction ([5]) includes additional terms in comparison 
with the useful approximation (see, for example, [18]) m 2 H ± = rn 2 A + m 2 ^ — e, where 
e = ?> G pmly / An 2 \/2 (m 2 / sin 2 (3 + m\j cos 2 (3) ln(Mg USY /m 2 op ) corresponds to the first 
term only in the expression for AA4, see ([5]). 



3 It is essential that (5) and (6) are valid for insignificantly different squark masses [TB]. The wave- 
function renormalization terms in AA4 [H [5] are omitted. For the case of large difference between the 
squark masses another decomposition of the effective potential must be used, see |17j . 
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Charged Higgs boson masses in the framework of typical MSSM scenario with strong 
mixing of the neutral CP-even/CP-odd states (known as the CPX scenario, see [TT] ) 
are shown in Fig.l. The two-loop corrections to A4 and the nonleading D-terms [I] are 
accounted for. In the CPX scenario under consideration the constraint on Msusy, ^t,b 
and /i is imposed as /i = 2 A tjb = 4M SUSY , everywhere in the following M SUSY = 500 
GeV. At moderate phases arg(/iv4 t 6 ) and small tg/3 the lightest neutral scalar mass is 
positively defined, see Fig. la, at relatively small starting from 20-30 GeV. The 
contours at the lightest neutral scalar mass m/ ll = 40 GeV, see Fig. lb, strengthen 
the limit on m H ± by approximately 20 GeV. At the same time at moderate values 
of the phase arg(/iv4 t fe ) of the order of 10-15 degrees the lightest neutral Higgs boson 
~ 40-50 GeV becomes poorly observable at LEP2 luminosities in the channel 
e+e" -> Zh 1 owing to the increase of the A scalar state CP-odd admixture in the hi 
state with the suppression of the ZZh\ coupling. The situation depends critically on the 
f^top value, at the top quark mass greater than 178 GeV and moderate tg/3 ~3-10 any 
value of could be possible down to values close to zero [19] . Small deviations from 
the CPX scenario with the A t ^ and fi parameters shifted by 100 GeV are illustrated in 
Fig. 2a and Fig.2b. In this case the region of smallest possible charged scalar masses is 
displaced to the large tg/3 ~ 40 — 45 interval. We are not considering in the following 
an extremely large tg/3 values exceeding 50 insofar as such values become to be poorly 
consistent with the Tevatron limits [20] on the t — > H + b decay channel. The t — > 
H + b decay probability is enhanced by tg/3 factor in the vertex and becomes poorly 
consistent with direct experimental data on the top quark pair production with the 
following decay t — > W + b. Direct LEP2 experimental data on the m#± production 
in the channels e + e~ — > cscs, cst~v t t + v t t^v t gives the limit m#± >56 GeV [2~T] . 
Note that the signal cross section is sensitive to potentially large one-loop corrections 
coming from superpartners and gluino [22] . Specific event charasteristics in the channel 
H ± — > W ± hi for a light charged Higgs produced in the top decay at the LHC and 
discovery possibilities were analysed in [23] . 

3. K°-K° mixing in the SM 

K°-K° mixing in the SM. The neutral kaon mixing in vacuum is characterized by 
the mass splitting of physical states Amis and by the parameter of CP violation e. 
In the framework of phenomenological approach to the mixing, see, for example, [24] . 
the effect is described by the nonhermitian Hamiltonian (mass matrix) of dimension 
2x2 and of the form My — ir^-, where the real M12 and the imaginary Ti2 parts 
of nondiagonal matrix elements define the value of Amis and the weights of CP- 
even/CP-odd states in the physical states of neutral mesons. The amount of CP 
violation in the phenomenological approach to the meson mixing is given by a small 
deviation from the ortogonality of the eigenvectors in the physical basis (given by 
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meson states without definite CP-parity which appear after the diagonalization of 2x2 
Hamiltonian) . The real and the imaginary parts of the off-diagonal matrix elements 
are defined by box diagrams (see Fig. 3) with the intermediate u, c, t quarks, and box 
diagram with the cut (c and t excluded), respectively. The SM mixing in the neutral 
kaon system appears as a consequence of mixing in the charged current sector defined 
by the CKM matrix [23], and a small value of the mixing is explained by the GIM 
mechanism [26J, see Fig. 3. 

The real part of the amplitude defines the mass splitting 
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taking into account the contributions of the order of {m t /mw) 2 and (m c /mw) 2 (the 
latter are small, so 1(6) ~ 1), and 1(6,65 £3) defines the contribution of "mixed" ct 
exchange diagrams (ct-boxes) [28, 29J: 
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(12) 



Here f K « 1.27 ^ 165 MeV is the decay constant, G F = 1.17 x 1(T 5 GeV" 2 , 
Vij are the CKM matrix elements, E>k ~ 1-0 is the nonperturbative QCD correction 
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("vaccuum insertion") and rji, i] 2 , rj 3 are factorized perturbative QCD corrections 
from gluonic "cross-exchanges" connecting inner and outer quark lines (see also [30] for 
two generations case). They are expressed by means of effective QCD couplings ratio 
taken at the corresponding quark mass scale in the power defined by the anomalous 
dimension. The leading logarithmic (LL) corrections are sensitive to fermion masses 
and the energy scale where we replace the IV-exchange amplitudes by the four-fermion 
ones, so corrections beyond LL are usually accounted for. For the kaon system most 
significant corrections are defined by rji because the main contribution is given by the 
cc box diagram. In the following we are using the NLO QCD corrections for rji, t] 2 , 
773 calculated in [31J. The NLO corrections to r] 2 and rj 3 contribute at the level of 
10% and are less important for the following consideration. We take rji = 1.3 using 
the parametrization of [31J at m c =1.3 GeV and A A j5=0.350 GeV, and r? 2 = 0.47, 
r/ 3 = 0.57. 

The contribution of the third generation quarks to e was obtained in [27J and [28] . 
The result ( TTUj) includes both the contributions of the cc-box and the combined ct-box; 
(jBJ) accounts for the tt-box and the combined ct-box yield to the neutral kaon mass 
difference. 



4. K°-K° mixing in MSSM II. 

In the absence of flavor- changing neutral scalar exchanges, two additional box diagrams 
with one or two charged Higgs boson exchanges appear in the MSSM II, see FigJH(.£/W 
diagram) and Fig 15] (HH diagram). In the calculation of HW exchange diagram using 
the t'Hooft-Feynman gauge, the unpysical scalar mode contribution was accounted for 
(Fig. 7). Using the low-energy approximation k 2 3> pf, where k/pi is the inner/(outer 
leg i) momentum, the kaon mixing parameters can be evaluated as follows: 

\ rn HW _ G F C H p K m K B K (, x 2 RoR , m \ m l p r> \ /-,o\ 

- [im. + m d )m H ReB 1 + — ReB 2 j , (13) 
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The notation Ch is used for the effective four-fermion coupling with scalar exchange, 
the GFermi analogue. Dimensionless functions Ihw-% Iiw-j, i — 1, 2, 3, 4; j = 1, 2, 
are the analogues of Inami-Lim-Vysotsky functions, appearing in the SM iy-exchange 
diagrams, and £j, i = 1, 6, are various mass ratios for particles of the internal lines, 

6 = (^) 2 , 6 = &) 2 , 6 = (^) 2 , e 4 = (^) 2 , £ 5 = (^) 2 , £e = (^) 2 - Symbolic 
expressions for Ihw-% an d Ihh-j and more details can be found in the Appendix. We 
are using r]\ = = rjj = 1.3, rj 2 = r] 5 = i] 8 = 0.47 and r] 3 = r] 6 = r/g = 0.57, 
defining the values of perturbative QCD corrections using [31]. Numerical contribtions 
of WW, HW and HH amplitudes to the neutral kaons mass splitting and to the e 
parameter are shown in Table 1 and Table 2. We define 

lF uot = 1 V™ + V™- 1 + VlT- 2 + Vl¥ riR , 
1 lLS ' " 2y/2W% w + W™- 1 + W™~ 2 + W™ 1 ' 

where V%", V™~\ V™~ 2 , V™ W? s w , W™~\ W™~\ W™ are the imag- 
inary and the real parts of various amplitudes multiplied by a corresponding effective 
factor 

VZ W = G 2 F -lmA, W% W = G 2 F -ReA, 



v hw-i = G F C H (m d + m 8 )m 2 H w hw-i = G F C H (m d + m s )m 2 H 



Vg^" 2 = m2 >< GFC » . l m Ba, = ^Ch^X, 

i5 8m^ sin 2/5 2 ' L5 8m^ sin 2/3 

We have already mentioned the importance of rji QCD correction factor in front of 
the cc box amplitude for the calculation of Am^ e. The NLO analysis of [31 J gives the 
uncerainty of rji about 20%: rji = 1.3 ± 0.2. We are using the central value T]i = 1.3. 
The values of ArriLs and e for the SM contribution at different values of r)\ can be 
found in Table 3. Relatively small HW and HH amplitudes have small influence on 
the full contrubution with the variation of rji factor. 



5. B° d -B° d and B° s -B° s mixings in MSSM II 

There is no conceptual difference in the analysis of the neutral B° s - and meson 
systems in comparison with the analysis of the K° system. Mixing of the states without 
definite CP-parity in vacuum with the mass splitting and CP violation defined by 
the CKM complex matrix elements and charged Higgs boson exchanges takes place. 
Numerical results can be obtained using formulae (jSJ), (TIP]) . ( TT3l) . (1T4"]) . ( |T5|) . (IIS]) and 
( TIT]) given above with the replacement of d or s index by the b index together with 
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redefinition of the decay constant, "vacuum insertion" and meson mass. Different quark 
masses and mixing parameter values lead to qualitatively different numerical results. 
The yield of SM cc-box amplitude to 5°-meson mass splitting Amf s is suppressed. 
Main contribution to mass splitting is given by the tt-box amplitude, the corresponding 
value of perturbative QCD correction factor t]b is not critically different from the r]i 
factor for the kaon system, rjs = 0.55 [32] . The u and c quarks provide important 
contributions to the diagrams with a cut, defining the imaginary part of nondiagonal 
Hamiltonian terms Ti2. Numerical values of Amg^ Amg s can be found in Table 4. One 
can observe that like for kaons, the contributions of HW and HH MSSM amplitudes 
to B° mass splittings are small. At m^± = 50 GeV the HW-box gives about 1% of 
B°-B° mass difference. The full splitting Am Bsd = Am^ + AmfJ^ + Amg^ is 
not significantly different from the SM result and is not in an excellent agreement with 
the experimental result. For example, in the SM we found B® mesons mass splitting 
Am Bd = 2.12 x 10~ 13 GeV and the experimental result is Am^J = (3.33 ± 0.03) x 
10~ 13 GeV. Additional contributions of HW and HH box amplitudes, see Table 5, 
parametrically dependent on the charged Higgs boson mass, can be used to restrict it. 

In order to evaluate the limits for m#± we used the approach described in [33J . Two 
statistical hypotheses are considered (1) charged Higgs boson exchanges contribute to 
the mass splitting (2) these exchanges do not contribute. The value of statistical un- 
certainty k is defined by both probabilities (a) the probability to reject the hypothesis 
(1) when it is true, (b) the probability to accept (1) when (2) is true. The relative 
uncertainty is defined as k = (a + b)/[2 — (a + &)], see details in [33], and converted 
to a number of standard deviations. For example, let us consider the case m#± =50 
GeV for I?°-mesons (see Table 4). Then the SM mass splitting is 2.11xl0~ 13 GeV 
and the negative contribution of charged Higgs exchanges is -0.017xl0~ 13 GeV. With 
the experimental accuracy 0.03xl0 -13 GeV indicated in the Table we get the relative 
uncertainty k (i.e. the probability of wrong decision about the existence of charged 
Higgs yield) of 32% which corresponds to 0.5<r level. If the experimental accuracy is 
taken to be 0.005xl0~ 13 GeV which is six times better, the relative uncertainty be- 
comes 2.5% corresponding to statistical significance of 2a (usually denominated as "the 
weak evidence" that the charged Higgs boson exchange contribution exists). With the 
experimental accuracy of the mass splitting measurement at the level of 0.003xl0~ 13 
GeV for 5^-mesons the restriction on the charged Higgs boson mass is m H ± > 1 12 
GeV with the statistical significance of 2a. The limit is stronger if the accuracy for 
5°-mesons is improved by one order of magnitude to the level of 0.01 xlO -13 GeV (see 
Table 4) giving m#± > 155 GeV at the 2a level. 

Note that for precise comparison of different amplitudes it is necessary to evaluate 
W ,H-boson exchange diagrams at small distances. Our approximation k 2 3> p 2 is not 
precise enough (see also [24]) from this point of view leading to an underestimate in 
comparison with the experimental mass splitting. More exact estimates can be found 



9 



in [5U ES] where results greater by a factor of 1.2 - 1.5 for the amplitudes have been 
obtained. 

The calculation of B° sd meson mixings is different from the case of i^°-mesons, 
because the phase of leading term in the imaginary part T 12 is approximately the same 
as the phase of the leading term of the real part M 12 , giving an additional suppression 
of the e parameter value by a factor {m c j m^) 2 . In this case the expression f TTTJ]) with 
the replacement V ts — > V t b cannot be used. 

6. Summary 

We consider the MSSM with explicit CP violation in the effective two-doublet 
potential and the Yukawa sector of type II (MSSM II) in the framework of the scenario 
opposite to the decoupling limit for heavy scalars. In the case of moderate SUSY scale 
(a few hundreds of GeV) relatively small Higgs boson masses and large Higgs bosons 
mass splittings of the order of masses the CP invariance of the effective potential is 
strongly violated. This scenario admits a relatively light charged Higgs boson m#± ~50 
GeV giving nonnegligible nonstandard contributions to the neutral mesons mixing. 

Using the low-energy approximation (k 2 3> p 2 ) we first evaluate the K° mass split- 
ting and the CP violation parameter given by (8), (13), (14), and (18). The yield of 
MSSM II amplitudes to K° mass splitting is extremely small far beyond the experi- 
mental precision and depends weakly on the possibly large tg/3=V2/v\ parameter, so 
the dominant contriution to Amis and mixing parameter e is defined by the standard 
VyVF-box amplitude. Analogously, for the B° s d - mesons the contribution of MSSM II 
B° — B° amplitudes is insignificant but much larger in comparison with the case of K° 
mesons system, especially for B° s mesons system at small charged Higgs boson mass. 
Charged scalar with mass 30-35 GeV is excluded if the precision of the B° mass splitting 
measurement equals 1%. Improvement of the experimental precision indicated in Table 
4 by about one order of magnitude excludes the charged Higgs boson mass m H ± < 112 
GeV for B^ mesons system and m#± < 155 GeV for B s mesons system. In principle 
the LEP2 limit could be improved by precise measurements on the 5-factories. 

Limitations on the charged scalar mass can restrict strongly the complex MSSM 
parameter space as illustrated by the two-dimensional contours in the variables (<p, tg/3) 
and (At,b, A*)j see Fig. 6a and 6b. Inside the contours shown in Fig. 6a which are 
generated for the CPX scenario (M SUSY =500 GeV) the lightest scalar mass m hl is 
positively defined. At = 60 GeV only the region inside the smallest contour in the 
vicinity of tg/5 =3 and zero phase is allowed, at the value m#± = 70 GeV this region 
slightly increases, and at m H ± = 80 GeV an additional allowed region appears at large 
tg/3 close to 40, which is broadening at m#± = 90 GeV (outer contours in Fig. 6a). 
At large tg/3 ~ 45 the lightest neutral scalar mass demonstrates high sensitivity 
to the values of (A tt b, /-0 parameters illustrated in Fig. 6b. The contours in Fig. 6b are 
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generated for =50 GeV and tg/3 =45, the inner contour at the value = 30 
GeV is increasing for charged scalar mass 50, 70 and 90 GeV (outer contour). Inside 
the contours exceeds 50 GeV. The allowed regions in Fig. 6b weakly depend on 
the particular choice of =50 GeV. So the light charged scalar with mass 50 - 60 
GeV in the scenarios with strong mixing of CP-even/CP-odd states and large tg/3 
does not leave much freedom for the choice of other complex MSSM parameters. 

We do not take into account the contributions of MSSM box diagrams with quark 
superpartners qi,R, chargino ^ and neutralino £°, defined by the trilinear MSSM La- 
grangian terms (for example, horisontal sides of a box are £°, and the 

exchange particles are Ql,r) ■ These Lagrangian terms include the 4x4 neutralino mix- 
ing matrix elements and 2x2 chargino mixing matrix elements, defined by the SUSY 
breaking mechanism, which are, generally speaking, complex numbers insofar as their 
phases cannot be removed by a redifinition of the spinor fields. The MSSM box contri- 
butions are dependent on a large parameter set involving sparticle masses and phases 
of matrix elements. Abundance of free parameters makes difficult the systematic esti- 
mates of MSSM contributions which can be very large or marginally small depending 
on a region of parameter space. This situation forces one to consider artificially simpli- 
fied MSSM scenarios even for a relatively simple observables in the states tt, W + W~ , 
see, for example, [36], which are accessible for direct reconstruction on the colliders. 
Explicit results for MSSM contributions to the meson mixing amplitudes in the class of 
models with radiative generation of fermion masses can be found in [37] . Following [37] 
the MSSM contributions increase strongly at large tg/3 around 40-50, specific choice 
of the phases and relations between the sparticles and gaugino masses. This scenario 
is not overlapping with our case when MSSM box amplitudes are expected to be not 
large at moderate tg/3 and degenerate sparticle masses. 

Our approximations are adequate for mixings at small external momenta pi in 
comparison with the momentum k in the loop. Evaluation of amplutudes at pi values 
which could be not small and probably give corrections of the order of 50% to Am^ 
and of the order of 10~ 2 to e ( see [27]) have been performed in [29], [31] for the general 
nonsupersymmetric two-doublet models with Yukawa sector of type I and II and also 
in [35] for the Yukawa sector of type III. Large loop momentum modifies Inami-Lim- 
Vysotsky functions [29]. The asymptotics of (21), (22) coincides with the asymptotics 
of modified functions in the limiting case m# — > corresponding to the SM WW 
box contribution only. 
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Appendix 



The simplest approximation for the one-loop box diagrams in the four-fermion limit 
m HW ^ ^ 2 ties U P the two quark-boson vertices to the one four-quark interaction 
vertex and the internal loop momentum k is removed from boson propagators. The 
lack of k powers in the denominators leads to a divergent amplitudes, so a cutoff is 
necessary at some momentum scale which is problematic to be meaningfully fixed. We 
are evaluating exactly the integrals over k which are convergent. The result is given by 
a dimensionless functions depending on the mass ratios for internal lines. In the SM 
they are known as Inami-Lim-Vysotsky functions which we denote in the case of HW 
and HH amplitudes by Iuw-i an d Ihh-j (see f lT5l) . ffT6]) and ffTTj) ). 

W TV-box, see Fig. 3 Details of calculation can be found in [27] . The integrands for 
various boxes are (using the t'Hooft-Feynman gauge and omitting a constant factor) 

d 4 k mttfnt A/ 1 

TP = C ' f W — (1Q) 

ww ~ l {k 2 + ml t ) 2 k 2 {k 2 + m 2 w ) 2 k e y } 

for cc- and tt-boxes with W exchanges, 

p d A k m 2 c m 2 m\y 1 



WW-2 - ^ 2 + m 2) (£2 + m ^)k 2 {k 2 + m 2 w ) 2 k G 

for ct-boxes with W exchanges, 



(20) 



d 4 k k 2 1 

jv a n k 

'ww-3 (k 2 + m 2 ) 2 (k 2 + m 2 w ) 2 ~ k 2 { 1 



for it-box with unphysical scalar, and 

d 4 k 1 
(k 2 + m 2 ) 2 (k 2 + m 2 w ) 2 ~ k 4 



for ii-box with W boson and unphysical scalar. 

The sum of these contributions gives ffTTj) and ffT2]) . In the SM the four-fermion 
approximation for PFW-box amplitude also converges, so keeping W propagators is 
not a principal improvement. 

HH-box, see Fig. 5 The integrand for cc and tt-boxes looks as 

p d 4 kWm 2 ct 1 

Ihh - 1 ~ (k 2 + ml t ) 2 (k 2 + m 2 H ) 2 ~ k 2 ' (23) 

Exact evaluation gives the dimensionless function 

T 6,2 + 6 26,26 



lHH-1 



(6,2 -£ 6 ) 2 (6,2-6 



In 6,5- (24) 
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For the ct-box the integrand has the form 

JP d 4 k k^k u m c m t J_ 

HH - 2 ~ (k 2 + m 2 c )(k 2 + m 2 )(k 2 + m 2 H ) 2 ~ k 2 ' 1 ] 

calculation of the integral leads to 

T - ft) + eM - g) + gg In 6 - gg In ft - gg jnft 

M " 2 (ft-ft) 2 (ft-ft) 2 (ft-ft) 

ftftft(ft-ft + 2ftlnft-2ftlnft) ^ 

(ft - ft) 2 (& - ft) 2 (ft - 6) ' 1 } 

where 6 = (^) 2 , 6 = (;^) 2 , 6 = & ft = ft = (^) 2 , 6 = (^) 2 - 

Unlike the SM PFVF-box the four-fermion amplitudes are divergent. 

HW -box, see Fig. 4 The evaluation is more or less simple in the t'Hooft-Feynman 
gauge when boxes with unphysical scalars (Fig. 7) must be accounted for. The inte- 
grands are 

d 4 k k^k v m\ t 1 
Irw - 1 = (k 2 + m\tf(k 2 + m\){k 2 + m 2 w )k 2 ~ F (2?) 
for cc- and tt-boxes, and 

t p d 4 k k^mcmtjmc + m t ) J_ , . 

w " 2 ~ k 2 (k 2 + m 2 c )(k 2 + m 2 ){k 2 + m 2 H )(k 2 + m 2 w ) ~ A; 4 1 j 

for ct-box. Dimensionless functions can be evaluated as follows: 

ft 2 2 ft(l - mft, 5 ) + ft 2 , 2 (lnft, 2 - 1) ~ ft, 2 ft 2 ~ 2ft, 2 ft In ft 

W_1 (l-e6)(6,2-e 6 ) 2 (6,2-l) 2 

, ft,2 + ft 2 (l + mft, 2 )-ft(l + lnft, 5 ) , , 

(l-£ 6 )(£l,2-£ 6 ) 2 (£l,2-l) 2 1 ^ 

and 

. ft 2 ft (In 6 - ft hi ft ) - Ci 2 Ce In ft + 6 In ft (6 + ft) + ft ft 2 ft In ft 

(1 - ftXft - ftXft - i)(ft - i)(ft - ftXft - 6) 

ftftg In ft - ft ft In ft - ft ft 2 In ft + gft In ft - ftg In 6 + 6ft In ft , v 
(l-ft)(ft-ft)(ft-l)(ft-l)(ft-ft)(ft-ft) ' 1 J 

The four-fermion amplitude is divergent. For the box diagram with unphysical scalar 
G in Fig. 7 

p d A kWm 2 ct 1 

Ih\V~3 = (ui i ™2 \ 2 /;„ 2 i _2 \V;„ 2 i ™2 \ ~ 72 (^1) 



(A; 2 + m2 t ) 2 (A; 2 + m 2 ^)(A; 2 + m 2 y ) A; 2 
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for cc- and tt-boxes, and 

j P d A k k^k v m c m t J_ 

w " 4 ~ {k 2 + m 2 c ){k 2 + m 2 ){k 2 + m 2 H ){k 2 + m 2 w ) ~ k 2 1 j 

for ct-boxes leading to the dimensionless functions 

T 6l,2 _ £ 1,2ft — Cl,2 — £l,2 l n ft,2 + 61,2^6 l n ft,5 + 61,2^6 

w ~ 3 " (ft-i)(ft, 2 -ft) 2 (ft, 2 -i) 2 

ft^ + 2ft, 2 ft In 6,2 ~ Ci, 2 eg ~ 2ft, 2 ft 2 In ft, 5 ~ ft 2 In ft 

(ft-l)(ft, 2 -ft) 2 (ft, 2 -l) 2 1 ^ 



and 



r g laftfe ~ 6) + ft 2 In 6 (6 - ft) + ftg(g In ft - In 6) 

W " 4 " (1 - ft) (ft - ft) (ft - l)(ft - l)(ft - ft) (ft - ft) 
ft 2 (ft 2 lnft - ft 2 In ft) + ft 2 ft 2 lnft(l - ft) - ftft(ft 2 ft hift + ft lnft) 

(l-ft)(ft-ft)(ft-l)(ft-l)(ft-ft)(ft-ft) ' 1 } 
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Figure 1: Charged Higgs boson mass in the complex MSSM II assuming that (a) the lightest 
neutral Higgs boson mass m/ ll is positively defined, (b) the lightest neutral Higgs boson mass 
TO/jj equals 40 GeV, as a function of tg/3 = 1)2/ v\. The sequence of contours is shown for the 
phases of two-Higgs-doublet potential effective parameters Xq, A 7 changing from zero (lower 
curve on both plots) to 180° (upper curve on both plots), with the step 10°. The CPX strong 
mixing scenario [11] is used at the scale Msusy =500 GeV. Below the contour for a definite 
phase value the lightest neutral Higgs boson mass (a) is not positively defined (b) is less than 



40 GeV. 
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40 42 44 46 48 50 40 42 44 46 48 50 

tg/3 

Figure 2: Charged Higgs boson mass in the complex MSSM II assuming that the lightest 
neutral Higgs boson mass m/ ll ~50 GeV at large values of tg/3 parameter. The curves are 
insignificantly sensitive to the limit imposed on m^ 1 . The sequence of contours is shown 
for the phases of two-Higgs-doublet potential effective parameters Xq, A 7 changing from zero 
(lower curve on both plots) to 180° (upper curve on both plots), with the step 10°. The 
strong mixing scenario close to CPX [11] is used at the scale Msusy =500 GeV (a) A t: b =890 
GeV, n =2000 GeV, (b) A t>h =890 GeV, fi =1900 GeV 
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Figure 3: The GIM mechanism for K° — > K mixing. 
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Figure 4: The i2W-box diagram for K° — > K° mixing. 
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Figure 5: The HH-hox diagram for K° — > K° mixing. 
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Figure 6: (a) Regions of lightest Higgs boson mass positively defined, (b) Regions 
of mh x >50 GeV, see the text for details. 
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Figure 7: The box diagram with unphysical scalar mode for K° — ► K° mixing. 
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Table 1: K° mesons mass spitting in the SM and MSSM II. The contributions of amplitudes 
shown in Figs. 3,4 and 5 to AM^s are indicated in the second, third and fourth lines. 



m,H±, GeV 


50 


100 


150 


200 


250 


300 


500 


Am^ s w x 10 15 , 


2.721 (Am^| = (3.482 ±0.013) x lO" 15 , ) 


Am^f x 10 19 , 


-14.22 


-4.44 


-2.19 


-1.31 


-0.87 


-0.63 


-0.24 


AmfJ" x 10 22 , 


61.50 


15.79 


7.15 


4.07 


2.63 


1.84 


0.67 



Table 2: Indirect CP violation for K° mesons in the SM (second line) and MSSM II 
(third line). 



ttih, GeV 


50 


100 


150 


200 250 300 500 


e 


ww x 10 3 


2.0523 (e^ = (2.232 ± 0.007) x 10" 3 ) 


£ 


: tot x 10 3 


2.0419 


2.0472 


2.0493 


2.0503 2.0509 2.0513 2.0519 



Table 3: Indirect CP violation and mass splitting for K° mesons in the SM at different 
values of rji. 





1.1 


1.2 


1.3 


1.4 


1.5 


Am e J x 10 15 , GeV" 2 


3.48 ±0.01 


sTs x 10 3 


2.23 ±0.01 


Am s L f x 10 15 , GeV" 2 


2.31 


2.51 


2.72 


2.93 


3.14 


e b L f x 10 3 


2.48 


2.25 


2.10 


1.88 


1.73 
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Table 4: B®- and i?!?-mesons mass splitting in the SM and MSSM II. Lines from two 
to seven indicate various contributions to AM T o 3 ■ 



inn, GeV 


50 


100 


150 


200 


250 


300 


500 


Am§7 x 10 13 , GeV 


2.11 (Am^ p = (3.33 ±0.03) x 10" 13 , ) 




Am^Y x 10 16 , GeV 


-17.05 


-11.99 


-8.70 


-6.57 


-5.14 


-4.13 


-2.09 


Amgf x 10 17 , GeV 


3.97 


2.87 


2.11 


1.59 


1.24 


0.99 


0.47 


Am™ w x 10 12 , GeV 


9.3 (Am e ^ = 1] 


lAt° H x lO" 12 , GeV) 


Am™ x 10 14 , GeV 


-8.21 


-5.67 


-4.08 


-3.07 


-2.39 


-1.92 


-0.97 


Amgf x 10 15 , GeV 


1.83 


1.31 


0.96 


0.72 


0.56 


0.45 


0.21 



Table 5: Relative contribtions to the mass splitting provided by HW-hoxes and SM 
WW-hoxes. R Bs d is the ratio of \Am™/Am^\. 



run, GeV 


50 


75 


100 


125 


150 


175 


200 


Am%Y x 10 13 , GeV 


2.11 (Am^f = (3.33 ±0.03) x lO" 13 , ) 


AmF x 10 16 , GeV 


-17.05 


-14.28 


-11.99 


-10.16 


-8.70 


-7.57 


-6.57 


R B d 


0.0081 


0.0068 


0.0057 


0.0048 


0.0041 


0.0036 


0.0031 


Am^ a w x 10 12 , GeV 


9.3 (Am^ p = 11.4+^ x 10" 12 , GeV) 


Ara^ x 10 14 , GeV 


-8.21 


-6.79 


-5.67 


-4.78 


-4.08 


-3.54 


-3.07 




0.0088 


0.0073 


0.0061 


0.0051 


0.0044 


0.0038 


0.0033 
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